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ABSTRACT 

The conditions that limit the possible excitation of ideal MHD axisymmet- 
ric ballooning modes in thin accretion disks are discussed. As shown earlier by 
Coppi & Coppi (2001a), these modes are well-localized in the vertical direction 
but have characteristic oscillatory and non-localized profiles in the radial direc- 
tion. A necessary condition for their excitation is that the magnetic energy be 
considerably lower than the thermal energy. Even when this is satisfied, there re- 
mains the problem of identifying the possible physical factors which can make the 
considered modes radially localized. The general solution of the normal mode 
equation describing the modes is given, showing that it is characterized by a 
discrete spectrum of eigensolutions. The growth rates are reduced and have a 
different scaling relative to those of the "long-cylinder" modes, commonly known 
as the Magneto Rotational Instability, that have been previously studied. 

Subject headings: MHD — accretion, accretion disks — magnetic fields — meth- 
ods: analytical 

1. Introduction 

The problem of identifying the processes which can produce significant transport of 
angular momentum outward in accretion disks has led to consider the plasma collective 
modes that can be excited in rotating plasmas where a magnetic field is present. In these 
plasmas the angular momentum is assumed to increase with the distance from the axis of 
rotation as in the case of Keplerian accretion disks. The most immediate approach to this 
problem is that of considering modes that are axisymmetric and are of the same type as 
those found originally for a long cylindrical plasma (Velikhov 1959; Chandrasekhar 1960; 
Balbus & Hawley 1991). 

Coppi & Coppi (2001a,b) pointed out that these axisymmetric modes when applied to a 
thin disk configuration become of the ballooning type (Coppi 1977) in the vertical direction, 
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i.e. they cannot be represented by a single Fourier harmonic, and acquire a characteristic 
oscillation in the radial direction. The radial wavelength is related to the distance over which 
the mode is localized vertically. Thus the problem of finding a localizing factor in the radial 
direction for the axisymmetric modes was left unresolved. As pointed out by Coppi & Coppi 
(2001c) the scaling of the growth rates of these ballooning modes is different from that of 
the original modes, which are appropriate for a long cylinder rather than for a thin disk, and 
are reduced by the condition of vertical localization. 

In the present paper we discuss the general solution of the two-dimensional equation that 
describes the axisymmetric modes in thin plasma disks where the magnetic energy density 
is much smaller than the thermal energy density. Two categories of modes are identified: 
the localized ballooning modes that have a discrete spectrum and the quasi-localized modes 
that have a continuous spectrum of growth rates and related radial wavelenths. These 
quasi-localized modes can be considered to be physically significant (i.e. to be sufficiently 
well- localized vertically) only for very large values of f3 l s = c 2 s /v 2 A where c s is the sound 
velocity and v A the Alfven velocity. As is well known, the ratio H/R is of the order of c s /v^ 
where = QR and Q = Q(R) is the local rotation frequency. 

When all the limitations are considered, even if a factor providing a radial localization 
for this kind of mode has not been identified yet, the rate of angular momentum transport 
that can be expected from its excitation does not appear to have the magnitude of those 
semi-empirical forms of the diffusion coefficients for the angular momentum transport that 
are used in current models of disks. 

In Section 2 the set of two-dimensional equations that give the spatial profile of axisym- 
metric modes is derived pointing out all the limits of the derivation. In Section 3 the vertical 
profile of the mode is shown to be described by a fourth order differential equation under 
realistic conditions. 

In Section 4, we give different analytical derivations of the localized solutions of this 
mode equation, corresponding to a discrete spectrum of growth rates and of radial wave 
numbers, in order to better describe the properties and the limitations of these solutions. 
The asymptotic decay of the modes in the vertical direction is proportional to a Gaussian. 

In Section 5 we introduce the class of "quasi-localized" modes whose asymptotic decay 
in the vertical direction is not proportional to a Gaussian but is algebraic and that have a 
continuous spectrum. 

In Section 6 the issue of the mode radial localization is discussed and the difficulty of 
constructing mode packets that could be radially localized is pointed out. We compare the 
ballooning modes to non-normal mode packets which may be formed from MRI modes to 
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achieve vertical localization. And finally the rate of angular momentum transport that can 
be expected from these modes is estimated and compared to that corresponding to currently 
adopted effective diffusion coefficients. 



2. Mode Profile Equation 



We adopt cylindrical coordinates and note that axisymmetric modes are represented by 
the perturbed toroidal velocity 

= Vt(R,z) exp(7 t) . 

In order to identify the properties of (R, z) we start by writing the total momentum 
conservation equation as 

A m = p ^ jU + v • Vv + v • Vvj + V (p + • B^ - ^B • VB ~ (I) 

where we neglect the R derivatives of all equilibrium quantities when compared to the R 
derivatives of the perturbation. Then, as is customary in the theory of ideal MHD modes, 
we consider the e^-Vx A m = component of Eq. (I). This is 

9 A 9 A n 

d~z AmR dR Amz ~ U 



and, specifically, 
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We note that, as shown by Coppi & Coppi (2001a,b), and the plasma displacement 
£ fl do not vanish when 7 /fi — > 0, while v R = 7 £k, v z = 7 £ 2 , and the compressibility 
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do. For the modes of interest the variation of (R, z) on z occurs on a scale distance that is 
considerably shorter than that for the variation of B z . This is consistent with the fact that, 
in the considered region, the B R component of the field is negligible. 

The frozen-in law can be written in the form 



^=(B-Vv + B-Vv)-B(V-v)-(v-VB + v VB) 



(3) 
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and the R component of it yields 



B R = B Z —£ R 
az 



9 u 9 u u 9 9 c 

By ~ B R ~ — By t R . 

dz z OR R z dRdz qR 



Thus Eq. (2) becomes 
Now, we write 



dz 2 
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-7o 2 P^V.^0. 

(4) 



considering the fact that for 7 = the frozen-in law (3) gives = — (dQ/dR)R^ R . Then 
the 4> component of Eq. (3) gives 

B, = B z ^-l-B,V-i 



and the component of Eq. (1) becomes 



P7o ( 7o |, + 2n&) ~ ±B,^ (b,^ -B^-i 



(5) 



Now the set of equations (4) and (5) needs to be completed by one that would relate 
V • £ to £ H and £ . Therefore we note that the adiabatic equation of state gives 

p ~ — pT (v ■ | 

and that by neglecting the z derivatives of the equilibrium quantities relative to those of the 
perturbations in the equation dA mz /dz = we obtain, for d£ z /dz = V • £ — d£ R /dR : 



<9z 2 







(6) 



with T = 5/3 and c 2 = Tp/p. 



Since, as will be shown in the following analysis, the excitation of axisymmetric modes is 
possible only if v 2 <C c 2 , we consider regimes for which this is the case and 7 2 <C c 2 \d 2 /dz 2 \. 
Then Eq. (6) reduces to 



dz 2 



(v.|) 
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c 2 dR q * c 2 B z dz^ 



(7) 
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This indicates that the terms involving V • £ can be neglected in Eqs. (4) and (5) and the 
set of equations that describe the axisymmetric modes reduces to 



p 7 o (Tel, + ^L) = ^ B zg^L (8) 



and 



d 2 r 



p ( 7o 2 |« - 2Q % i, + 2QQ'RQ - —B 2 Z 



( 
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i R + 




0. (9) 



dz 2 
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In order to assess the influence of the modes that can be excited on the effective rate 
of transport of angular moment in thin accretion disks, a number of factors have to be 
considered. Thus the most evident factors to be taken into account in looking for the 
possible solutions of Eqs. (8) and (9) are: 

i) The mode growth rate 

ii) The vertical and radial characteristic scale distances 

iii) The value of the threshold for the mode onset 

iv) Whether or not a mode is contained (localized) vertically within the disk and 
whether appropriate non-normal modes (packets) can be constructed 

v) The physical features that may localize the mode radially. 

In particular, it is well known that the rate of transport which can be produced by excited 
plasma modes can be assessed in a rudimentary way from an effective diffusion coefficient 
that includes both the mode growth rate and its characteristic wavelenths. Therefore modes 
with very short wavelengths may be less effective that modes with longer wavelengths even 
if the growth rates of the latter modes are smaller. 

For the sake of simplicity, we do not consider at first the condition that a mode should 
be vertically localized within a disk nor the criterion that it should have the lowest threshold 
fl c possible for its onset (Q > Q c ). In particular, we refer to modes that are propagating in 
the vertical direction, have vertical wavelenths X z = 2n/k z that are much smaller than the 
height of the disk, 2H, and involve the central value of the particle density p = p(z = 0). 
We assume also that the modes are propagating radially. Therefore 



(R, z) 



exp [ik z z + ik R (R - R )] 
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where R = R is the radius at which all equilibrium quantities entering Eqs. (8) and (9) are 
evaluated, and k z H 3> 1. Thus Eq. (8) becomes 

(7o 2 + ^L)l* = -7o2fi|« 

where J\ z = k 2 z v 2 Az , v\ z = B 2 /(Att Po ), and Eq. (9) 

( 7o 2 + c4) k 2 l R - 3il 2 k% - 2n % i,k 2 z = 

where k 2 = k\ + /c 2 and we have taken 2Q'QR = — 3f2 2 as appropriate for a Keplerian disk. 
The resulting dispersion relation is 

To 4 + To 2 (2*L + fi 2 |) = (s^ 2 | - *l) • (io) 

We see that the instability threshold is 3f2 2 = k 2 v Az , implying very high values of f3 s = (? s /v 2 Az 
as k 2 > k 2 ~^> 1/H 2 , and that the maximum growth rate (7 relatively close to Q) is found 
for k\ <^ k 2 . This is the limit in which the Velikhov (a.k.a. MRI) instability is found. The 
question then arises whether the assumptions that have to be made to arrive at Eq. (10) can 
make these types of modes good candidates to produce the needed rate of transport in thin 
accretion disks. Therefore, in the next section we look for the modes that can be contained 
in the disk vertically and have the lowest threshold possible. 



3. Contained Modes 

Now we analyze the characteristics of the modes that are localized near the center of 
the disk over a distance smaller than H where the density p{z) can be approximated by 

z 2 \ 



and H ~ H . 



These modes, as was demonstrated by Coppi & Coppi (2001a,b), are characterized by 
being evanescent in the vertical direction and oscillatory in the radial direction and are 
represented by 

£ R (R, z) ~ l R {z) exp [-ik R (R - i?„)] 

where £ R (z) is an even or odd function of z such that £ fl (|-2| — > H) — > 0. Two other important 
points are that, contrary to the case where Velikhov modes with the highest growth rates 
can be found, 



k 2 R > 



1 d 2 ~ 



i R dz 2 
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and the relevant growth rates (Coppi & Coppi 2001c) scale as 

7c 2 ~ |f < & • (12) 

Therefore axisymmetric modes that are contained within the disk have a characteristic 
growth rate that depends on the scale height of the disk and is well below that of the 
MRI instability. 

We note that, as shown by Eq. (7), |V • £| <C \d£ a /dR\. Therefore ik R £ a ~ —d^ z /dz, 

z 

Uz) ~ ik R [ Uz')dz' , (13) 



and £ z (z) has the opposite parity of £ R (z). Moreover, since £ z (z) has to be localized, like 
£ R (z), within the disk the solution £ R (z) is subject to the condition 

3 

Uz')dz' = . (14) 



The problem to be solved is simplified if we note that Eq. (8) reduces to 

d 2 ~ 



(15) 



by observing that 7 2 ~ v 2 Az /H 2 and that the modes we look for are localized within — H < 

z < H : 



H 



1 d 2 ~ 



Under the same conditions Eq. (9) reduces to 



> 1 . 



dz 2 



3Q 2 z 2 ~ 

(k 2 R v 2 Az - 3Q 2 + 7o 2 ) G + - 2 7o^ " v 2 Az — Lr 



if 

Az Tz r 



Thus the equation for the mode amplitude is reduced to one of the fourth order, that is 



dz 2 



2 2 
R U AZ 



3n 2 + 7o 2 ) i R + 



3tt 2 z 2 ~ 



d 2 z 



m G v ****** 



7X + ^]i,-- 



Referring to Eq. (15) we see that £ ^ only when 7o ^ 0, and in this case the condition 
that — > if) — > implies that 



(16) 



— oo — oo 
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Therefore when 7 7^ 0, £ R (z) has the two constraints (14) and (16) in addition to that of 
vertical containment £ R (\z\ — > H) — > 0. 



Now we define 



with z = z/A z and obtain 



, 2 _ 3Q 2 




Az 



dz 2 



d 2 ~ ~ (k 2 7 2 



A; 2 4\ 7 2 # 2 ~ 



It is evident that the solution is localized over a smaller distance than H if k 2 H 2 ^> 1. 
We note also that k\jk 2 — 1 ~ l/(k H ) ~ v Az /c s <C 1, and we consider the asymptotic limit 
(12). Therefore we introduce two additional dimensionless quantities involving the radial 
wave number and the growth rate 

A,-(§-lW 



% 2 



*o 2 " 3 



For k R ~ A; + A ~ —25k R H . The implied dimensionless form of the normal mode 
equation is thus 

dz 2 [dz 2 



r2 -z 2 +A)U-z) 



+ i%m = o . 



(17) 



Finally we note that the adopted linearized approximation is valid to the extent that 
\£ g \<A z , implying \£ RM \y/k H < y/WJk~ , that is 



6 



< 



VazH 



k 



where \£ RM \ 
A < k H . 



max|^ fl |. We note also that the condition k 2 jk\ — \ <C 1 implies that 



4. Solutions of the Normal Mode Equation 



The simplest case to consider is that of marginal stability that is 7 2 
(17) reduces to 

~ d 2 1 ~ 

— -z 2 + A = 



0. Then Eq. 
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with the boundary condition that £ R (z) vanishes at large z. The solutions of this equation 
are given by 

~ _ ( z 2 

in{z) = £ R , H n (z) exp I - — 

for integer n > 0, where H n are the Hermite polynomials, with the corresponding eigenvalues 
A = 2n + 1. 

We note that in order to comply with the condition that £ z be a localized function of 
z, represented by Eq. (14), only odd-parity £ R (z) solutions are acceptable. Thus the lowest 
eigenfunction is represented by 



£r(z) = 6^ exp 



-z 2 



and 



L(z) = (-i (k H ) 1/2 | fl ,o) exp 



The corresponding eigenvalue is then A = 3 and k 2 R = k 2 [1 — 3/(k H )]. 

In order to analyze a case where the solution extends beyond the values of z 2 /H 2 for 
which the approximation (11) for the particle density is valid, we consider the complete 
profile p = p exp(—z 2 /H 2 ). Then the equation for the marginally-stable modes is 



exp 



2 \ 



H 2 ) 



d 2 ~ 







(18) 



k 2 

exp(-l 2 ) - -f 



and we see that k R cannot be too small relative to k in order to have the needed turning 
points within the main body of the disk. Moreover, if the solution is to be well-localized 
within the disk the condition (k H ) 2 ^> 1 corresponding to (3 S ^> 1 is necessary. An example 
of this is illustrated in Fig. 1 where k H = 16. To illustrate this it may be useful to rewrite 
Eq. (18) as 

f » (l)+ (db^ f » (i)=0 ' 

where z = z/H . The solution given in Fig. 1 has been obtained numerically and compared 
to the one that would be found analytically replacing exp(— z 2 ) by 1 — z 2 for an equal value 
of k H . 

Next we consider the case of unstable (7 ^ 0) modes that can be found for A > 3. We 
note that the relevant solutions of Eq. (17) are subject to the conditions that ^ in addition 
to t; z and £ H vanish at large z. Therefore it is convenient to recast the equation in terms of 
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the variable Y oc considering that Y' oc £ 2 and Y" oc £ K given Eqs. (13) and (15). Then 
we rewrite Eq. (17) as a set of two second-order differential equations 



dz 2 



d 2 z 



dz 2 



Ir + if - A) In = %Y 



(19) 
(20) 



and note that this set of equations corresponds to the following matricial equation 



On 
22 



Y 



A 7o 
7o 



is. 

Y 



where the operators On and O22 are easily identifiable from Eqs. (19) and (20). This has a 
discrete spectrum of the eigenvalue matrix, when the conditions that both £ R and Y are of 
a given parity and that (, R (z 2 — > 00) — > and Y{z 2 — > 00) 



0. 



Then we refer to the equation for Y 



that lends itself to derive the following quadratic form for £ z 



(21) 



i" 



S2 



A 



+ 7 2 



6 



(22) 



This is obtained multiplying Eq. (21) by d 2 Y/dz 2 and integrating over (—00, +00). 

The localized solutions for z large are given, as in the case where % — > 0, by 

d 4 ~ 9 rf 2 ~ 
—r-rY — z ——rY ~ . 

Thus we have six constraints for the solution of Eq. (21): normalization, two parity con- 
ditions, and three unwanted asymptotic limits, for large z, out of the expected four. Six 
constraints lead to two eigenvalue conditions: only specific values for both A and 70 are 
allowed, as anticipated on the basis of Eq. (22). 

The exact solutions that we have identified for Eq. (21) confirm the general conditions 
and characteristics that we have outlined. In particular the lowest eigensolution is 



(23) 



Y = exp ( -— ) oc ^ 
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and corresponds to 

z 2 \ dY 



and 



L = Ik,o (1 - exp (~y) > 

as shown in Figs. 2 and 3. We note that, for Y(z) given by Eq. (23), Y"" = (z 2 - 1)Y" + 
2Y + AzY' and find A = 5, % 2 = 2 from the condition 3 - 7 2 - A - (5 - A)z 2 = 0. Thus the 
relevant growth rate is 

and 



«*» = (25) 



The next eigenfunction Y = z exp(—z 2 /2) corresponds to £ z = £z, (l — z 2 ) exp(— z 2 /2) 
and £ fl = £ H ,oz(3 — z 2 ) exp(— z 2 /2). The relevant set of eigenvalues is A = 7, 7o = 6. 
Moreover, there are two even-parity solutions at A = 9, two odd solutions at A = 11, three 
at 13 and 15, and so on up to higher and higher eigenvalues, as indicated in Table 1. 

We observe that the eigenfunctions with the highest values of 7 2 corresponding to the 
same value of A (that is, the most unstable solutions) are, in fact, given by 

Y(z) = H n e~*<* (26) 

A = 5 + 2n (27) 
7 2 = (n+l)(n + 2). (28) 

The problem can be dealt with more simply by noting that the Fourier transform of the 
relevant solutions satisfy the following second-order equation 

ti^fa) + Ah ^k f ^ + ^ - ^ + 2 - *£) f = 

whose solutions corresponding to Eq. (26) are 

Y(k,) = k 1 iexp (-% 



For these solutions we have 



2 + 3n + n 2 - 7 2 + (A - 5 - 2n) k\ = 
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which gives the conditions (27) and (28). Combining the two, we obtain the dispersion 
relation 

7 2 = i (A-3)(A-l) (29) 
where A > 3 takes integer odd values. 



5. Quasi-localized Modes 

In addition to the well-localized modes with the discrete spectrum discussed in the 
previous section, there is a class of modes to consider which decay algebraically rather than 
exponentially as a function of z and that have a continuous spectrum (see Figs. 4 and 5). 
Therefore we may classify them as quasi-localized modes and, in practice, consider only 
those which have relatively fast power-law decays such as Y ~ l/z a with a > 5 as physically 
relevant. 

The asymptotic solutions for \z\ — > oo of the mode equation that depend algebraically 
on z are given by 

d 2 ~ 2 ~ 

d¥ Y -^-°> 

and are of the form 

Y oc 2 V*±ViM (30) 

The "practical" condition a > 5 corresponds to 7^ > 30 which can be considered to exceed 
the asymptotic limit (12) for which the mode equation (21) has been derived, unless k 2 H 2 
has extremely large values such that 30 < 7^ <C k H . 

The four asymptotic solutions of the fourth-order equation are: decaying exponential, 
decaying power-law, growing exponential, and growing power-law. The last two are of course 
unacceptable. Excluding only the growing solutions, we are left with five conditions on a 
fourth-order equation. This leads to one eigenvalue condition, corresponding to 70 being a 
continuous function of A. 

Since the discrete solutions are a special case of the continuous spectrum formed by 
a linear combination of Gaussian-decaying and power-law decaying asymptotic forms, we 
consider that Eq. (29) remains valid for continuously varying values of A. Then we note that 
7 2 > 30 corresponds to A > 13. 

We may observe that when A and 7o are related by Eq. (29), Eq. (21) can be factored 
as follows: 

( d 2 d A - 1\ / d 2 ~ d ~ A - 3 ~\ , x 

\^-z— + }[—Y + z—Y + Y =0 31 

\dz 2 dz 2 J \dz 2 dz 2 J K ' 
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This factorization separates the growing from the decaying asymptotic solutions. Thus 

^*«p (-y) +C2 i(A^ 

The discrete spectrum corresponds to values of A for which c 2 = 0. A relevant case is 
illustrated in Fig. 4. 

We note that the second-order equation resulting from Eq. (31) can be solved analytically 
in general. Introducing the variable u = —z 2 /2, we have 

4 ? C" + G- M )^M-(^) ? M = ° 

which is the standard differential equation for the confluent hypergeometric function iiq, 
also known as the Kummer function: 

v f-\ v f A ~ 3 1 **\ 
Y even {z) - I — - — , -, -— J 

Y odd (z) = z X F X 7,1 

For discrete choices of A and particular parities, these Kummer functions simplify to the 
localized solutions of Eq. (26), but in general they have power-law decays at large z, and 
solutions exist of either parity at the same eigenvalues. 



6. Observations 

For the modes that have been presented, it is clear that A will have to be limited by 
the condition that 7„ ~ v 2 Az /H 2 <C Q 2 . Therefore unless v 2 Az /c 2 is exceedingly small very few 

1/2 

modes can be fitted in the interval 5 < A <C k H ~ (3 S . Thus a problem which arises is 
that of identifying a physical factor which can localize radially the normal modes considered; 
a mode packet, constructed by varying A over a significant interval, is difficult to envision. 

We note also that, by comparison with the "long cylinder" modes, the "cost" of localizing 
a mode vertically between z = —H and z = +H and creating the needed turning points 
in the appropriate equation involves (i) the requirement that k R be significant and related 
to the scale of vertical localization A 2 ~ (H /k o y/ 2 with k R A z 3> 1 and k R ~ k ; (ii) a 
considerable reduction of the growth rate 7 relative to Q and its dependence on the scale 
distance H . If a physical factor providing the radial localization of the discrete modes given 
by Eq. (21) can be found it is reasonable to expect that the resulting growth rate will depend 
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on a second scale distance and that it should be depressed further relative to the expression 
(24). In order to resolve the issue of the radial localization of the mode, the best option may 
be that of attempting to superpose modes centered around a range of values of R , taking 
into account that the corresponding values of k R , j , and A z vary as a function of R . 

The argument could be made that by taking very small wavelengths, requiring that 
very large values of f3 s be considered, the MRI dispersion relation (10) could be used with 
k 2 R <C k 2 (see Fig. 6). Then "large" growth rates would be obtained, but wave packets would 
have to be constructed to ensure that the resulting non-normal mode would be localized 
vertically. However, we can argue easily that a packet can spread rapidly, losing its initial 
localization property. In particular we consider an interval of values of k z around k z = k° z 
where the growth rate given by Eq. (10) is maximum. Therefore in this interval, 



lo (k z ) ~ 7o (£;°) + 



l d 2 7o 
2 dkl 



(k z k° z ) 2 , 



(32) 



where d 2 ^ /dk 2 < 0, and d%/dk z = for k z = k° z . Then we construct the packet 



and we obtain 



J dk z exp |-A 2 (A; 2 - k° z ) 2 + \^t(k z - K) 2 + i(k z - k° z )z^ 
-y (K)t + ik° z z + ik R (R - R ) 



(33) 



cxp 



\l/ oc exp 



2 dkl 1 



+ 7 t + ik° z z 



(34) 



where A 2 H 2 . It is evident that the packet begins to spread for k° z (<9 2 7 / 'dk 2 )t ~ 2A 2 k° z . 

In order to estimate the rate of angular momentum transport that these packets may 
produce we derive the flux of angular momentum from the quasilinear theory of the uncon- 
fined normal modes that are the mode packet components. Clearly, 



T.j = R 



Considering the perturbations from the equilibrium state, the flux Tj can be separated 
into an average and a fluctuating part: Tj = (Tj) + f j where (f j) = 0. In particular 
(Tj) = J2 k (Tjk) and we note that the particle flux associated with these modes (r p ) is null 
as (pv R ) = given the fact that V • v = and p — 0. Then 



(Tj k ) = p (% k v* Rk ) - - (B, k B* Rk ) + c.c. 



(35) 
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where v Rk = % k Z, Rk , = -Q'R£ Rk + % k ^ k , B Rk = i(k z B z )£ Rk , and B^ k = i(k z B z )^ k . For 
G k = (T Jk /(pR)) we obtain 



G k — 27 ofc 





to 




; 



-fl'R + 2Q 



ilk 



and the relevant effective diffusion coefficient can be defined as 



1 - 



2^ - ll k 



SVR u 2 Az + 7o 2 fc 



(36) 



(37) 



We note that the transport of angular momentum, produced by modes of the considered 
type for which V • v = and evaluated from the relevant quasilinear theory assuming the 
adiabatic equation of state, is not accompanied by an increase of thermal energy. This is 
an important feature that differentiates the effects of these modes from those of a viscous 
diffusion coefficient. 

Using the dispersion relation (10) we can verify that T> c r is positive. To assess the order 
of magnitude of £> e ff we may argue, as is usually done in plasma physics, that, at saturation, 
\^ R \ ~ l/k° z ~ l/k . Thus we have 



V, 



cff 



10 
U2 



n 



H 

% 



(38) 



Since the existence of these modes requires that ft > 1 we may conclude that this kind of 
mode packet should not produce a rate of transport of angular momentum that is comparable 
with the one represented by the Shakura-Sunyaev coefficient (Shakura & Sunyaev 1973) 



V ss ~ a s 



,c,H 



(39) 



This is frequently used with ad hoc values of the numerical coefficient a ss that exceed 10 2 . 

It is evident that the mode packets have the intrinsic limitations of not being the normal 
modes of the system when considering their possible excitation. In particular we may consider 
their growth to be limited to the case where 



A - 



2 dkl al 



where a H is a reasonable numerical coefficient such as 3. A numerical derivation of the mode 
packet (34) deriving d 2 ^ /dk 2 z from the dispersion relation (10) is given in Fig. 7. We note 
that the packet spreads rather slowly as a function of time, as A cff = [A 2 — (d 2 ^ / dk 2 z )t} 1 / 2 . 
On the other hand the initial spread A has to be relatively broad as the dispersion relation 
(10) limits the range of values of k z for which Eq. (32) is valid, as indicated in Fig. 6. The 
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consequence of this is that the packet spread can reach the region where the Alfven velocity 
varies significantly relative to its central value even at the outset. 

When dealing with normal modes that are contained within the disk, given the discrete 
spectrum of both k R and 7 that characterize them and that the vertical eigenf unctions are 
of the ballooning type, the standard quasilinear theory cannot be applied to arrive at an 
estimate of T> c r. We may in fact use the expression (38) which can be arrived at by standard 
qualitative considerations, and conclude that the relevant values of D c r would fall well below 
those estimated from Eq. (39). On the other hand, the greatest difficulty that the ballooning 
modes have is that on envisioning a process by which they can be localized radially. 

It is a pleasure to thank P.S. Coppi, with whom the present work was started, for his 
continuing interest, and L.E. Sugiyama for comments. This work was sponsored in part by 
the U.S. Department of Energy. 



A. Higher Branch Solutions 

We note that additional branches of solutions to Eq. (21) also exist at higher A at exact 
intervals of 4, but with lower corresponding values of 7„ . For completeness, the solutions for 
branch m > are given by: 

m ( A — 3 1z 2 \ 

Y even (z) = a m,i( A ) iF 1 li-m+ — — , -, - — j 

V> , s v^, n / A-l 3 z 2 \ 
Y odd {z) = z o m>i (A) i-ti \ i — m-\ — , -, -— 1 

i=o ^ ' 

where the coefficient functions are given recursively by: 

... (i-m-l)(A + 4i-4m-7) ... 
a '" AA) = «(A + «-4 m -2) a -- l(A) 

= (i- m -l)(A + 4i-4m-5) 
, v ; i(A + 2i-4m-2) ' v ; 

with the base case given by normalization: a mj o(A) = 6 mj o(A) = 1. The dispersion relation 
for the higher branches is 

_ 2 _ (A -Am- 3) (A - 4m - 1) 
To - ^ 

with A > 3 + 4m. 
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Fig. 1. — Corrections to the marginally stable solution. The solid line is the analytic 
marginally-stable solution for the approximation (11) of the density profile, giving A = 3. 
The overlapping dotted line is a numerical solution for the mode in the case of the density 
distribution p = p exp(—z 2 /H 2 ) that gives A ~ 2.9. Here \fk Q H = 4, a relatively low value, 
is chosen to emphasize the difference between the two solutions. 
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^R ^R 

Fig. 2. — Vertical profiles of the marginally-stable mode (left) and of the lowest unstable 
mode (right) for y/k H = 7. Here k = y/3Cl/v Az , and H is the characteristic scale distance 
for the density vertical profile near z = 0. In the figure z is scaled to H . The equilibrium 
density profile is shown as a grayscale background. The value of k H should be sufficiently 
large that the mode is well-localized over < \z\ < H , e.g. £(z = H )/£(z = 0) 0.1. 
Note that the mode is localized over the width A z ~ H /y/k H and that (fc -ff ) 1//2 ~ fi x J A \ 
Therefore the corresponding values of f3 s are very large. 
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Fig. 3. — Two-dimensional streamline patterns of the plasma displacement for the 
marginally-stable mode (left) and the lowest unstable mode (right). Solid versus dotted 
lines denote the sign of rotation. Length units in the figure are scaled to H . As noted in the 
text, the radial variation of the mode is faster than its vertical variation, which is in turn 
stronger than the vertical variation of the density profile. Here y/k H = 7. 
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Fig. 4. — Profiles of a localized even mode (upper) and of an odd quasi-localized mode 
(lower) for A = 13 and ^ = 30. The decay of the even mode is proportional to a Gaussian 
while that of the odd mode is proportional to z~ 5 , where z = z/H. Here y/k H = 7. 
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Fig. 5. — Graphic representation of the dispersion relation for the axisymmetric ballooning 
modes. The marginally-stable modes (70 = 0) and the unstable modes (70 > 0) corresponding 
to discrete and odd values of A are indicated by dots. Additional branches of solutions with 
similar behavior exist at multiples of 4 at higher A. Clearly the highest branch corresponding 
to the highest growth rates for a given value of A is the meaningful one. The continuous 
spectrum of quasi-localized modes is indicated by a dotted line, but in practice they cannot 
be considered as physically acceptable unless 7o > 30, a value that exceeds the asymptotic 
conditions under which the mode equation (21) has been derived. 
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Fig. 6. — Dispersion relation of the MRI modes. This plot shows the growth rate of the MRI 
modes as a function of k z when k R = 0. The maximum growth rate of 7 ~ 0.850 occurs at 
k z ~ 0.6k , and the region of instability is bounded by k 2 R + k\ < k%. Here k = y3Q/v Az . 
The short dashed line gives the approximation (32) to the dispersion relation. The long 
dashed line is a sample wavepacket with A = lOA;^ 1 centered around the peak growth rate, 
shown for comparison with the accuracy of the approximation. Clearly, narrower spectra for 
k z could be chosen, but they would correspond to larger values of A. 
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Fig. 7. — Spreading of non-normal mode packets. An initially contained mode packet formed 
from a range of k z centered around the peak growth rate of the MRI will spread over time 
due to the dispersion relation. The top panel gives an initial packet with a width A = 10k~ l 
at t — 0, and the bottom panel shows the same packet at j Q t = 25. Here k = y/3fi,/v Az , and 
the width A is chosen by the criterion that the mode packet undergoes minimal deformation 
as a function of time while having the maximum growth rate. The length is given in units 
of /c" 1 , so one can see that the mode packet is unacceptably spread even at t = for 
k H < 30 (corresponding to (3 S ~ 1000), and even higher values of j3 s are necessary to retain 
containment through many orbit times. 
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Table 1. Normal Mode Solutions 



A 


To 2 


Y{z) 
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2 




7 
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ze^ 2 ' 2 
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12 


(1 - z 2 ) e- 22 ' 2 
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(l + lz 2 )e- s2 / 2 


11 


20 


(z - |z 3 ) e-* 2 ' 2 


11 
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(z + 2z 3 ) e~~ z2 l 2 


13 
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(l-2z 2 + \z A )e~~ z2 / 2 


13 
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(l + z 2 - \z A ) e-* 2 / 2 


13 
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(l + ±z 2 + ±z*)e-* 2 / 2 
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